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Abstract 

A general procedure for studying finite-N effects in quantum phase transitions of finite systems is pre- 
sented and applied to the critical-point dynamics of nuclei undergoing a shape-phase transition of second-order 
(continuous), and of first- order with an arbitrary barrier. 

An important issue concerning quantum phase transitions in mesoscopic systems is to understand the 
modifications at criticality due to their finite number of constituents. In the present contribution we study 
this question in connection with nuclei exemplifying a finite system undergoing a shape-phase transition. We 
employ the interacting boson model (IBM) ^Jj which describes low-lying quadrupole collective states in nuclei in 
terms of a system of N monopole (s) and quadrupole (d) bosons representing valence nucleon pairs. The model 
^) is based on a U{6) spectrum generating algebra and its three dynamical symmetry limits: C/(5), SU{3), and 
CnI 0(6), describe the dynamics of stable nuclear shapes: spherical, axially-deformed, and 7-unstable deformed. A 
Q I geometric visualization of the model is obtained by an intrinsic energy surface defined by the expectation value 
D " of the Hamiltonian in the coherent (intrinsic) state [2j [3] 
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where = (1 p'^)-^^^[f3 cos7dJ -I- /3sin7 {d!^ + d^_2)/\/2 + s^']. For the general IBM Hamiltonian with one- 
^ t and two-body interactions, the energy surface takes the form 

£;(/3,7) = Eq+N{N (3^)-^ [a(3'^ -h(3^ 005^,-1 + c(3'^] . (2) 

■ The coefficients Eo,a,b,c involve particular linear combinations of the Hamiltonian's parameters [4]. The 
quadrupole shape parameters in the intrinsic state characterize the associated equilibrium shape. Phase tran- 
sitions for finite N can be studied by an IBM Hamiltonian involving terms from different dynamical symmetry 
chains [3] . Several works have followed this route in numerical studies of finite-N effects at criticality [5-8] . In 
1-^^ ' the present contribution we consider an (approximate) analytic-oriented approach to this problem [9-11]. 
I , The nature of the phase transition is governed by the topology of the corresponding intrinsic energy surface 

^ ' which serves as a Landau's potential. In a second-order phase transition, the energy surface is 7-independent 
^ , and has a single minimum which changes continuously from a spherical to a deformed 7-unstable phase. At the 
^ critical-point a = b = and the energy-surface acquires a flat behaviour (~ /?"*) for small /3 (justifying the use of 
J> , a square- well potential in the E(5) critical-point model [12 ). This is the situation encountered in the C/(5)-0(6) 
' phase transition where the critical Hamiltonian is given by 

' 1 r 1 r 1 

; H^r, = e7id + ^^[rf^-dt"(s^)^J [rf-rf-s^J , e={N-l)A. (3) 

Here d^ = {—l)^d-^ and the dot implies a scalar product. Hcri is 0(5)-invariant and involves the U{5) term 
(the d-boson number operator), and the 0(6)-pairing term. The intrinsic energy surface of Hcri has the form 

E{(3) = ^AN{N -1) + AN{N -l)f3*{l+ (3^)-^ . (4) 

As shown in Fig. (la). The global minimum at /3 = is not well- localized and E{f3) exhibits considerable 
instability in (3. Under such circumstances fluctuations in /3 are large and play a signiflcant role in the dynamics. 
Some of their effect can by can be taken into account by means of variation after projection of states of good 
0(5) D 0(3) symmetry tLM from the coherent state in Eq. ([T]) 
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These r-projected states form the ground band (^ = 1) and interpolate between the C/(5) ground state, \ N,nd = 
T = L = 0) = \s^) at /3 = and the 0(6) ground band, \ N,a = N,t, L) at /3 = 1. The matrix element of 
Hcri © in these states defines the r-projected energy surface which can be evaluated in closed form [9] 
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Figure 1: Energy surfaces of the critical IBM Hamiltonian Hcri © with iV = 5 and A ^ 1. (a) Intrinsic 
energy surface E{(3), Eq. ^ [solid line], and its approximation by a square- well potential [dashed line], (b) 0(5) 
projected energy surface i?^=i,r=o(/3), Eq. ^ whose global minimum is at = 0.314. 
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2) have approximate wave functions of the form 



jj-^ ^YiQ states (O, and are given by S\"^' = 
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with Pp = ■ S — /?^(s^)^] and J\fp a known normalization. Explicit expressions for quadrupole rates involving 
r-projected states can be derived. For example, with the E2 operator T{E2) = rftg + s^d, 



B{E2- C = 1; r + 1, i' = 2r + 2 ^ ^ = 1, r, L = 2r) 
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Similar expressions for the excited-band energies E^^l^{l3) and interband E2 rates are available As seen 
from Fig. (1), in contrast to E{f3), the lowest 0(5) projected energy surface E^=i^r=o{P) supports a well-defined 
global minimum at a certain value of (3. As shown in Table 1, using this effective /3-deformation in the r-projected 
states provides accurate analytic estimates to the exact finite-N calculations of the critical IBM Hamiltonian 
which in-turn capture the essential features of the E(5) critical-point model [12] present in ^■^'^Ba. 



Table 1: Excitation energies (in units of £'(2+^ — 1) and B(E2) values (in units of B{E2; 2+^ 0+q) — 1) 
for the E(5) model [12], for several N=5 calculations and for the experimental data of ^■^^Ba. The finite-N 
calculations involve the exact diagonalization of the critical IBM Hamiltonian (Hcri) [Eq. ([3|)], r-projected 
states, L+ , [Eq. © with /J^ = 0.314], U{5) and 0(6) limits of the IBM. 
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Figure 2: Energy surfaces of the critical IBM Hamihonian, Eq. ([9|), with n = 0.2 and N = 10. (a) Intrinsic 
energy surface £'(/?) = -E(/3,7 = 0), Eq. (fTO| . [sohd hne]. The unmixed L = and L ~ 2 levels are shown for 
illustration, (b) E{P) [solid line] as in (a), unmixed L = Q [dashed line] and L = 2 [long dashed line] projected 
energy surfaces, El{(3) = e\^\p), Eq. and the lowest L — eigenpotential [dot-dashed hne], E''j~}oiP), 

whose global minimum is at = 0.591. 



In a first-order phase transition the intrinsic energy surface has two coexisting minima which become de- 
generate at the critical-point. This is the case in the U{5)-SU{3) phase transition where the critical Hamiltonian 
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involves the U{5) pairing and the SU{3) quadrupole terms. The associated intrinsic energy surface 
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has two degenerate minima, at /? = and at (/3 — -j^tI ^0)- shown in Fig. (2), the barrier separating 
the spherical and prolate-deformed minima is extremely small and the resulting surface, E{P) = E(P,j = 0), 
is rather flat. This behaviour motivated the use of a square-well potential in the X(5) model [13 . Particularly 
relevant are states of good 0(3) symmetry L projected from the intrinsic state [/3, 7 = 0; N) of Eq. ([1]), 
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where Pl{x) is a Legendre polynomial with L even and r^''(/3) is a normalization factor. The i-projected states 
[/3;iV, L,M) interpolate between the C/(5) spherical ground state, |s^), with = t = L = 0, at /3 = 0, and 
the SU{2)) deformed ground band with (A, /i) = (2A^, 0), at /3 = v^- The matrix element of the Hamiltonian in 
these states define an L-projected energy surface which can be evaluated in closed form [10] 
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Here s[^j^ , ^2^1!^ ^2^1! respectively, the expectation values of fis 



s's, {s^)^s , hsfi-d in \I3\N,L,M) and 



are given by ^j^^ = r^^,(/3)/r^^(/3), with S)^^^ = ^I^l ^l!!"'' and Y.'^'l = {N - 1)S['2 " ^2!!' ■ As shown in 
Fig. (2b), the L = projected energy surface, E^j^q{(3), no- longer exhibits the double minima structure observed 
in the (unprojected) intrinsic energy surface. Instead, there is a minimum at /3 > 0, a maximum at /? = 0, and a 
saddle point at /3 < 0. The L — 2 projected energy surface, E^j'^J^{j3)^ has a minimum at a larger value of /3 > 0, 
and a flat shoulder near /3 = 0. The different behaviour of E^j^J^{l3) and E^j^Jq{(3) can be attributed to the fact 
that the L = 2 state is well above the barrier and hence experiences essentially a flat-bottomed potential. In 
contrast, the two minima in the intrinsic energy surface support two coexisting spherical and deformed L — 
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Table 2: Excitation energies (in units of E{2j) = 1) and B(E2) values (in units of B{E2; 2^ 0^) = 1) for 
the X(5) critical-point model [13], for several N=10 calculations, and for the experimental data of ^'''^Sm. The 
finite-N calculations involve the exact diagonalization of the critical IBM Hamiltonian [Eq. ([9])], L-projected 
states [Eqs. (HH) with f3 = 0.591], U{5) and SU{3) limits of the IBM. 
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states which are subject to considerable mixing. This mixing can be studied by means of a 2 x 2 potential energy 
matrix, Kij , calculated in the following orthonormal L = states 

l*i> = , \^2)^{l~rl^)-"^[\l3,N,L = Q)-r,^\s^)) , 

K.J - (*.|i?e,.|^,) , n2^{s''\f3;N,L = Q) = [mr\^=°\/3)r'^' . (13) 



The derived eigenvalues of the matrix serve as eigenpotentials, E'j^Q{P), and the corresponding eigenvectors, 
l^^^^g), are identified with the ground- and first-excited L = states. The deformed states \(3]N,L,M) of 
Eq. (jlip with L > are identified with excited members of the ground-band with energies given by the L- 
projected energy surface E^j^\(3), Eq. (12). As shown in Fig. (2b), the lowest eigenpotential E^j~]q{/3) has a 
global minimum at a certain /3 > 0. Using this value in the L-projected states leads, as seen in Table 2, to 
faithful estimates to the exact finite-N calculations of the critical IBM Hamiltonian (notably for yrast states), 
which in-turn capture the essential features of the ^(5) critical-point structure [13] relevant to -'^^^Sm. 

In a general first-order phase transition with an arbitrary barrier, the intrinsic energy energy surface of 
Eq. ([2]) satisfies a,b > and — Aac at the critical-point. For 7 = it can be transcribed in the form [11] 

E,„iP) = Eo+cN{N-l)fiP) 

f{/3) = 13^1 + (3^)-^l3 - Pof . (14) 

As shown in Fig. 3, Ecri{P) exhibits degenerate spherical and deformed minima, at /3 = and (3 — — > Q. 
The value of [3q determines the position (/? = /?+) and height {h) of the barrier separating the two minima in a 
manner given in the caption. To construct a critical Hamiltonian with such an energy surface, it is convenient 
to resolve it into intrinsic and collective parts [^ ; 

Hen = Hint + He ■ (15) 

The intrinsic part {Hint) is defined to have the equilibrium condensate | /3 = /3o, 7 = 0; TV), Eq. ([1]), as an exact 
zero-energy eigenstate and to have an energy surface as in Eq. (14) (with Eq = 0). Hint has the form 



Hint — ho 



/3ostdt + ^7/2 (dt^t)(2) , /3oJs+ y772(Jd)(2) , (16) 



and by construction has the L-projected states \(3 = /3q; N, L) of Eq. pip as solvable deformed eigenstates with 
energy E — Q. It has also solvable spherical eigenstates: | A^, Ud = t = L = 0) = \s^) and | A^, = r = 3, L = 3) 
with energy E — and E — 3/i2 [/3o(^ ~ 3) -I- 5] respectively. For large A^ the spectrum of Hint is harmonic, 
involving quadrupole vibrations about the spherical minimum with frequency e, and both /3 and 7 vibrations 
about the deformed minimum with frequencies e = = PqN ^ e-y = 9{1 + Pq)~^ C/j, where the last inequality 
holds in the acceptable range < Po < l-i. All these features are present in the exact spectrum of Hint shown 
in Fig. 4, which displays a zero-energy deformed (K = 0) ground band, degenerate with a spherical (rid = 0) 
ground state. The remaining states are either predominantly spherical, or deformed states arranged in several 
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Figure 3: The IBM energy surface, Eq. at the critical-point of a first-order phase transition. The position 
and height of the barrier are /3+ — and h — /(/3+) — j (^-'^ + a/I + /Sq ) respectively. 

excited K = bands below the 7 band. The coexistence of spherical and deformed states is evident in the 
right portion of Fig. 4, which shows the Ud decomposition of wave functions of selected eigenstates of Hint- The 
"deformed" states show a broad nj, distribution typical of a deformed rotor structure. The "spherical" states 
show the characteristic dominance of single Ud components that one would expect for a spherical vibrator. 

The collective part {He) of the full critical Hamiltonian, Eq. (15), is composed of kinetic terms which do 
not affect the shape of the intrinsic energy surface. It can be transcribed in the form [?| 



H, = C3 C, 



-0(3) 



Grid + C5 C, 



-^0(5) 



471. 



En 
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where N = fid + fis- Here Cq denotes the quadratic Casimir operator of the group G as defined in [J^j . Table 3 
shows the effect of different rotational terms in He- For the high-barrier case considered here, (/3o = 1.3, 
h = 0.1), the calculated spectrum resembles a rigid-rotor {E ^ ai\fL{L + 1)) for the cs-term, a rotor with 
centrifugal stretching {E ^ umL^L -|- 1) — hN[L{L + 1)]-^) for the C5-term, and a X(5)-like spectrum for the 
cg-term. In all cases the B(E2) values are close to the rigid-rotor Alaga values. This behaviour is different 
from that encountered when the barrier is low, e.g., for the U{b)-SU{?>) case discussed above, corresponding to 
/3o — 1/2-^/2 and h « 10""^, where both the spectrum and E2 transitions are similar to the X(5) predictions. 



Table 3: Excitation energies and B(E2) values (in units as in Table 2) for selected terms in the critical Hamil- 
tonian, Eqs. ■ . (fTTl) . The exact values are for N = 10, /?o = 1-3 and cl//i2 coefficients indicated in 
the Table. The entries in square brackets [...] are estimates based on the L-projected states, Eq. pTjl . with 
(3 = 1.327, 1.318, 1.294, determined by the global minimum of the respective lowest eigenvalue of the potential 
matrix, Eqs. ([T8| . (fT9|l . The rigid-rotor and X(5) [13] predictions are shown for comparison. 
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Figure 4: Left potion: spectrum of Hint, Eq. (fT6| . with /12 = 0.1, /3o = 1.3 and N = 10. Right portion: the 
number of d bosons (rid) probabihty distribution for selected eigenstates of Hint- 



Considerable insight of the underlying structure at the critical-point is gained by examining the 2x2 
potential energy matrix, Kij, of Eq. (|13p . which reads 
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K22 = So + (1 - ri^y [Eillil3) + 2ca'N{N - l)r^ 
Apart from a constant shift, e')'^\(3) is the L-projected energy surface of Hcri (fTS]) . and is given by 
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Here ^2^'' , d[^j^ , D^j^ and 5*2^^ denote the expectation values of nsfid, fid, fid{nd — 1) and hg {hg — 1) respectively 
in 1/3; N, L, M) (fTTj) . As seen in Table 3, by determining the value of (3 in the L-projected states from the global 
minimum of the lowest eigenvalue of the potential matrix (jlSp . one obtains accurate finite-N estimates to the 
energies and E2 rates at the critical-point. 

Part of the work reported was done in collaboration with J.N. Ginocchio (LANL). This work was supported 
by the Israel Science Foundation. 
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